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Abstract 

Motivated by the problem of reducing the peak to average power ratio (PAPR) of transmitted signals, we 
consider a design of complementary set matrices whose column sequences satisfy a correlation constraint. The 
design algorithm recursively builds a collection of 2 t+1 mutually orthogonal (MO) complementary set matrices 
starting from a companion pair of sequences. We relate correlation properties of column sequences to that of the 
companion pair and illustrate how to select an appropriate companion pair to ensure that a given column correlation 
constraint is satisfied. For t = 0, companion pair properties directly determine matrix column correlation properties. 
For t > 1, reducing correlation merits of the companion pair may lead to improved column correlation properties. 
However, further decrease of the maximum out-off-phase aperiodic autocorrelation of column sequences is not 
possible once the companion pair correlation merit is less than a threshold determined by t. We also reveal a design 
of the companion pair which leads to complementary set matrices with Golay column sequences. Exhaustive search 
for companion pairs satisfying a column correlation constraint is infeasible for medium and long sequences. We 
instead search for two shorter length sequences by minimizing a cost function in terms of their autocorrelation and 
crosscorrelation merits. Furthermore, an improved cost function which helps in reducing the maximum out-off-phase 
column correlation is derived based on the properties of the companion pair. By exploiting the well-known Welch 
bound, sufficient conditions for the existence of companion pairs which satisfy a set of column correlation constraints 
are also given. 
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I. Introduction 

Complementary sequence sets have been introduced by Golay [1], [2], as a pair of binary sequences with the 
property that the sum of their aperiodic autocorrelation functions (ACF) is zero everywhere except at zero shift. 
Tseng and Liu [3] generalized these ideas to sets of binary sequences of size larger than two. Sivaswamy [4] and 
Frank [5] investigated the multiphase (polyphase) complementary sequence sets with constant amplitude sequence 
elements. Gavish and Lemple considered ternary complementary pairs over the alphabet {1, 0, —1} [6]. The synthesis 
of multilevel complementary sequences is described in [7]. These generalizations of a binary alphabet lead to new 
construction methods for complementary sets having a larger family of lengths and cardinalities. However, all these 
studies focus either on the set complementarity or on the design of orthogonal families of complementary sets. 
Correlation properties of column sequences of the complementary set matrix (i.e., the matrix whose row sequences 
form a complementary set) have not been considered. 

In [8]— [10], a technique for the multicarrier direct-sequence code-division multiple access (MC-DS-CDMA) 
system [11], [12] that employs complementary sets as spreading sequences has been investigated. Each user 
assigns different sequences from a complementary set to his subcarriers. By assigning mutually orthogonal 
(MO) complementary sets to different users, both multiple access interference and multipath interference can be 
significantly suppressed. Similar to conventional multicarrier systems, one of the major impediments to deploying 
such systems is high peak-to-average power ratio (PAPR). We have stressed in [13] that correlation properties of 
column sequences of complementary set matrices play an important role in the reduction of PAPR. In this work, 
we search for ways of constructing complementary set matrices whose column sequences satisfy a correlation 
constraint. 

For orthogonal frequency-division multiplexing (OFDM) signaling, Tellambura [14] derived a general upper 
bound on the signal peak envelope power (PEP) in terms of the aperiodic ACF of the sequence whose elements 
are assigned across all signal carriers. He has shown that sequences with small aperiodic autocorrelation values can 
reduce the PAPR of the OFDM signal. By generalizing earlier work of Boyd [15], Popovic [16] has demonstrated 
that PAPR corresponding to any binary Golay sequence (i.e., a sequence having a Golay complementary pair) is 
at most two. This has motivated Davis and Jedwab to explicitly determine a large class of Golay sequences as a 
solution to the signal envelope problem [17]. Here, we consider sequence sets which are characterized by both their 
complementarity and a desired column correlation constraint. 

We describe a construction algorithm for the design of 2* +1 MO complementary set matrices of size 2 t m by 



2 t+p+1 in Section III, where t and p can be any non-negative integer, and m is an even number. The construction 
process is based on a set of sequence/matrix operations, starting from a two column matrix formed by a companion 
sequence pair. These operations preserve the alphabet (up to the sign) of the companion pair. In Section IV, we 
illustrate how, by selecting an appropriate companion pair we can ensure that column sequences of the constructed 
complementary set matrix satisfy a correlation constraint. For t = 0, companion pair properties directly determine 
matrix column correlation properties. For t > 1, reducing correlation merits of the companion pair may lead 
to improved column correlation properties. However, further decrease of the maximum out-off-phase aperiodic 
autocorrelation of column sequences is not possible once the correlation of the companion pair is less than a 
threshold determined by t. We also present a method for constructing the companion pair which leads to the 
complementary set matrix with Golay column sequences. 

In Section V, an algorithm for searching for companion pairs over length m sequences of a desired alphabet is 
described. However, exhaustive search is infeasible for medium and long sequences. We instead suggest finding 
companion pairs with a small, if not minimum, column correlation constraint. In Section VI, by exploiting 
properties of the companion pair, we convert the problem into a search for two sequences of length m/2 with 
low autocorrelation and crosscorrelation merits, a long standing problem in literature (e.g. see [18]— [21]). We 
further derive an improved cost function and show how it leads to reduced achievable maximum out-off-phase 
column correlation constraint. Sufficient conditions for the existence of companion pairs which satisfy various 
column correlation constraints are also derived. We conclude in Section VII. 



Throughout this paper, sequences are denoted by boldface lowercase letters (e.g., x), their elements by 
corresponding lowercase letters with subscripts (xq), boldface uppercase letters denote matrices (X), and calligraphic 
letters denote either sets of numbers or sets of sequences {X). 

A. Correlation functions 

Let a = (ao,ai, ...,a n _i) denote a sequence of length n with G C, < i < n — 1, where C is the set of 
complex numbers. The aperiodic and periodic ACFs of a are 



II. Definitions and Preliminaries 
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TABLE I 
Correlation Merits 
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where a* denotes the complex conjugate of a, and © n denotes modulo-n addition. It follows that 

P a (l) = A a (l) + A a (n-l), 0<Z<n-l. (3) 

Let b = (60, 61, b n -\), where 6« G C, < z < n — 1. The aperiodic and periodic crosscorrelation functions of 
a and b are defined, respectively, as, 
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Table I lists the correlation function parameters which are commonly used to judge the merits of a sequence 
design, and are termed correlation merits. X a > "^a > an( ^ &a ^ autocorrelation merits, and A^ b , A^ b , S£ b , and 
S^ b are common crosscorrelation merits. For example, it is well known that binary m-sequences satisfy = 1, 
and the sequences with \£ < 1 are called Barker sequences. Furthermore, a small value of can significantly 
reduce the PAPR of OFDM signals, if the elements of a are assigned across all carriers [13], [14]. 

B. Complementary sets and column correlation constraints 

A set of m sequences {ai,a2, ...,a m }, each of length n, is called a complementary set if 

m 

5]4(l) = 0, l</<n-l. (6) 
i=i 

When m = 2, the binary sequence set {ai, a2} is called a Golay complementary pair and a«, i = 1,2, are Golay 
sequences. They are known to exist for all lengths n = 2 Q T0 /3 26 7 , where a,(3,j > [22], [23]. Complementary 
sets {bi, b2, b m } and {ai, a2, a m } are mates if 

TO 

J> ai)bi (0=0, (7) 
i=l 



for every I. The MO complementary set is a collection of complementary sets in which any two are mates. Let n 
denote the MO complementary set consisting of k complementary sets each having m complementary sequences 
of length n. For binary sequences, k cannot exceed m [23], that is, the maximum number of mutually orthogonal 
complementary sets is equal to the number of complementary sequences in a set. Hence, M™ ra is called a complete 
complementary code of order m [24]. 

A complementary set can be represented using a complementary set matrix 
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where T denotes the matrix transpose. Row sequences, = (cj j0 , c^i, Cj^, ■ ■ ■ , Cj )Tl _i), 1 < i < m, are 
complementary sequences, that is, YLlLi ^r 4 (0 = 0, 1 < / < n — 1. The main focus of this paper are properties of 
column sequences, Cj = (cij, C2 j- , c^j , • • • , c mj ), < j < n — 1. 

A set of k mutually orthogonal complementary sets, {Ci,C 2 , ...C^}, form a MO complementary set matrix 
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and its column sequences are denoted as Uj, < i < kn — 1. An upper bound on an autocorrelation merit of column 
sequences is termed the column correlation constraint. If there exists at least one MO complementary set matrix 
satisfying a given column correlation constraint, then this constraint is called an achievable column correlation 
constraint. For example, let 

A„ = max {Au., < i < kn - l] , (10) 

if A„ < X A , then n is called a MO complementary set matrix satisfying a column correlation constraint X A , and 
X A is an achievable column correlation constraint. We also consider column sequences which are Golay sequences, 



i.e., Golay column sequences. 



C. Companion pair 

Let a be a sequence of length m, where m is an even number. We define a sequence b as a companion of a if 



is a complementary set matrix which consists of m/2 complementary pairs. C is called a companion matrix and 
(a, b) is called a companion pair. 



TABLE II 

Sequence Operations 
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D. Operations and extensions 

1) Sequence and matrix operations: In Table II, we list the following sequence operations: reversal, negation, 
complex conjugation, concatenation, interleaving, and inner product. Furthermore, we introduce two sequence 
reshaping functions /$(■) and f c (-) defined on sequences of even length. 

Let C = [cij] and D be two matrices of equal dimensions, then C* = [c* •] and — C = [— Cij]. C D is the 
matrix whose ith row sequence is obtained by interleaving ith row sequences of C and D. CD denotes the matrix 
whose ith row sequence is the concatenation of ith row sequences of C and D. 

(v) (v) 

Let us also define a sequence set 1Z C ' , which is a collection of row sequences of matrix tv c ' recursively 



constructed from a sequence c, as follows, 
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Let K%d = n c' u n d ) which consists of 2 V+2 sequences of length 2 v m. For example, 7Z C = {c, — c}, 
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= {d, -d}, and = {c, d, -c, -d}. 



2) Complementary set matrix extension operations: We describe two operations for extending complementary 

set matrices, namely, length-extension and size-extension. 

Lemma 2.1 [23]: Let {a, b} be a complementary pair, then {b*, —a*} is its mate, and both {ab*,b(— a*)} and 

{a (g) b*, b (8) (—a*)} are complementary pairs. 

Proof: See Appendix A. ■ 
Lemma 2.2: Let {ai,bi}, {a2,b2j-, {a m ,b m } be m complementary pairs of length n. Then, 

{ai, bi, a2, b2, ...,a m ,b m } is a complementary set of 2m complementary sequences. 

Proof: E£i(^a i (0+^b i (0)=0, l</<n-l. ■ 
Lemmas 2.1-2.2 imply that, if a complementary set consists of m/2 complementary pairs, the sequence length 

can be recursively doubled as follows. Let 
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be, respectively, an m by ra( p ) complementary set matrix and its mate, where m is an even number, and 
{r^, r^}, {r^\ rjj^}, ...{r^i_ 1} r$?} are assumed to be complementary pairs. A complementary set matrix C( p+1 ) 
of dimension m by 

„(p+i) = 2n (p) can be constructed recursively as either, 

c (p+i) = c (p) D (p) (15) 
c (p+i) = c^fg-D^. (16) 



or 



We term dl5t and (11 6b length-extension operations. 



Lemma 2.3 [3]: A MO complementary set matrix M|^ 2n can be constructed recursively as either, 
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Proof: Refer to the proof of Theorem 12-13 in [3]. 
We term ( 1171) and ( 1181) size-extension operations. 
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III. Construction of complementary set matrices from a companion pair 
A. Recursive construction 

Let X[m) denote a sequence set which consists of all length m sequences whose elements are from the alphabet 
X . We summarize a recursive construction of a MO complementary set matrix M 2 * 2t (p) with elements from X, 
where n& = 2P +1 , m is an even number, and t,p = 0, 1, 2, .... 

Step i: The construction starts from a companion pair Co and ci which are in X{m). They form an m by 2 
companion matrix 
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Step 2: By employing the length-extension operation p times, we extend C^ ) to an m by = 2 P+1 complementary 
set matrix C( p \ and its mate D^ p ) constructed from Eq. dT4l form a MO complementary set matrix 
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Step 3: Starting with (f0 , we can construct the MO complementary set matrix M 2 ** 2 « (j>) by repeating the 

size-extension operation t times, where p, t = 0, 1, 2, .... 

In this paper, we will alternately use either "the constructed MO complementary set matrix" or, simply, 
M^tl^ 2 t (*>) when referring to the above constructed MO complementary set matrix. 



B. Companion pair design and properties 

Proposition 3.1: Let us arrange the elements of Co = (ci,o, C2,o> •••c m ,o) into m/2 arbitrary pairs, e.g., (c Xj o, c^o)- 
Then, its companion sequence ci = (01,1,02,1, ...c m) i) can be constructed as either 
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Proof: Let us assume ci is constructed 
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Here, rl 0) = (c xfi ,c yfi ) and = 
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— c* ), respectively, the xth and the yth row sequence of form a 


complementary pair, since 
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Based on Lemma 2.2, C^°> is a complementary set matrix consisting of m/2 complementary pairs. Hence, c and 
Ci form a companion pair. ■ 
Example 3.1: /*(cq) is a companion of cq, since 
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is a companion matrix. 
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The companion pair has the following properties, 

Property 1 (Commutative property): If Co is a companion of ci, then ci is also a companion of Co- 
Proof: In Eq. dTTt. C is still a companion matrix when the column sequences a and b are switched. ■ 
Property 2 (Inner product property): If Co and ci form a companion pair, then Co • ci = 0. 
Proof: Based on Eqs. dm and d22l . J2^=i c i,o c i,i = 0- ■ 
Corollary 3.1: Binary sequences Co and ci form a companion pair, if and only if cq • ci = 0. 
Proof: From Property 2, Co • ci = is for any companion pair. Furthermore, if Co and ci are binary sequences 
such that Co • ci = 0, there must exist m/2 pairs of (x, y) satisfying c x fiC x ,\ + Cy,oCy,x = 0, where 1 < x < m and 
1 < y < m. Hence, row sequences of C^ ) can be arranged into m/2 pairs, where each pair ri°^ = (c x ,o, c x ,i) and 



,(o) 



(cy,o,Cy,i) satisfies Am (I) + A to) (I) = 0, 1 < I < 2. Consequently, C^ ) is a companion matrix. 



C. Column sequence properties 



Lemma 3.1: All column sequences of the constructed complementary set matrix M^ n(p) = [C^ D^ p ^] are in 
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{±cq,±ci}, where (cq,ci) is the companion pair for the construction. 



Proof: Let 
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For C^ p+1 ) constructed by the length-extension d 1 5I ». we have, 

c„ (p) < i < nW - 1 
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It follows that column sequences of C( p+1 ) are equal to, or are a negation of, column sequences of C^ p \ Since 
column sequences of C^ ) are c and ci, all column sequences of C^ v \ p = 0, 1, 2, are in TZ^J yCl - In addition, 
for C( p+1 ) constructed using the length-extension (fTol l. the interleaving of the corresponding row sequences doesn't 
change the column sequences and, thus, its column sequences are also in TZ^ jCl . Any column sequence of D( p ) 
can be found in C^ p+1 ^ and, consequently, it is in TZ^ Cl as well. ■ 

Lemma 3.2: All column sequences of the constructed MO complementary set matrix 2 *n» WQ m T^i]^, 

where (co, ci) is the companion pair and t,p = 0, 1, 2, .... 

Proof: Let {uf\o < i < r = 2 2m r>) } denote column sequences of 2 t n<p) • ^ ne size-extension (fT7b 

implies, 
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where < i < r. Based on Lemma 3.1, uf\ i.e., column sequences of (p) = [C^ D^'] , are in 7^ )Cl - 



From (O, O. and (EU, we have that 



uf G^£ ClJ 4 = 0,1,2,... (29) 
When M|t^ 2 *r!» is constructed using size-extension (fT8b . the proof is analogous. ■ 

IV. Properties of the constructed complementary set matrix 

In this section, column correlation properties of the constructed MO complementary set matrix are related to 
ACFs of the companion pair. We illustrate how to satisfy a column correlation constraint by selecting an appropriate 
companion pair. We also construct the companion pair which leads to complementary set matrices with Golay 
column sequences. Since number of zeros in an energy-normalized sequence can affect its PAPR (see e.g. [6], [9]), 
we also discuss the number of zeros in column sequences at the end of this section. 

A. Column correlation properties 

Theorem 4.1: MO complementary set matrix M|* * 2t w satisfies a column correlation constraint, if and only if 
the companion pair (co,ci) is selected so that all sequences in TZ^ )Cl satisfy the constraint. 

Proof: The proof is a direct consequence of Lemma 3.2. ■ 
Corollary 4.1: Complementary set matrices and constructed from a companion pair (co,ci), satisfy a 
column correlation constraint, if and only if cq and ci satisfy the constraint. 



Proof: The proof follows by setting t = in Theorem 4.1. ■ 
The minimum achievable column correlation constraint for the constructed MO complementary set matrix 
^2*m 2*n(p) * s a function of its size and alphabet and can be expressed as follows, 

^min(t, m) = min |max | A c : c G lZc}, Cl j : (co, cx) is a companion pair and cq, ci G A'(m) j , (30) 

where A is any autocorrelation merit. The following Lemma 4.1 is a. key in relating column correlation constraints 
to the correlation of the companion pair. In particular, it leads to the minimum achievable column correlation 
constraint X^ in (t,m). 

Lemma 4.1: The ACF of any column sequence u- 4 ^ of M|«t* 2* («o can ^ e ex P resse( i i n terms of the ACF of Co 
or ci , recursively, as follows, 

2A w {l)±A w {a-l) i 0<l<s; 
A^ + i ) (l) = < • ' (3D 

±A iv)(l -s), s<l < 2s. 

where '+' holds for j = i + r or j = i + 3r, '— ' holds for j = i or j = i + 2r, when size-extension H71 is used; 
'+' holds for j = 2% + 1, '— ' holds for j = 2i, when size-extension dT8l is employed; u. ^ G {co,cx, — Co, — ci}; 
< i < r = 2 2v+1 n(P\ s = 2 v m, and v = 0, 1, 2, t - 1. 

Proof: Eq. i3ll can be derived based on (128b . Note that, u 4 - ^ G 7£co!ci = { c o> c i) — c (b — c i} an d the negation 
of a sequence doesn't change its ACF. ■ 

Similar recursive equations can be found for periodic ACFs based on Eq. ©. 

Proposition 4.1 (A sufficient condition for S A ): Let S A < S A , S a < S A , and A Co (0) = A Cl (0) = E, where 
(co, cx) is a companion pair. Then, a sufficient condition for M 2 ,^ 2 < w to sat i s fy the column correlation constraint 
S A is 

S^>A t s£ + 2 t ~ 1 (2 t -l)E. (32) 

Proof: Let {u| , < i < 2 2 * +1 n^} be the column sequences of M 2 ,^ 21li(p) . Clearly, G {c , cx, — c , — cx}. 
Then, based on d3Tl . we have that 

S A = max{s A t) ,0<i<2 2t+1 n^} 

{2'm-l "| 
l^ u w(0l 5 °< * < 2 2m n^ I 

< 4*S^ + 2*- 1 (2*-l)E, (33) 
for t = 0, 1, 2, .... Hence, d32l is sufficient for 5 4 A > S„ which proves the proposition. ■ 



Proposition 4.2 (A sufficient condition for X A ): Let X A o < X A , A A < X A , and ^4 Co (0) = A Cl (0) = E, where 
(cq, Ci) is a companion pair. Then, a sufficient condition for M 2 ^ 2t (p) to satisfy the column correlation constraint 



X A is 



X A > max {(2* - 1)E, (2 t+1 - 1)A^} . (34) 

Proof: Let {uf\o < i < 2 2t+1 n^ pS> } be the column sequences of M 2 *^ 1 2 t n w implies 

A„ = max{A^ t) ,0<i <2 2W >>} 

= max||A u ( t )(Z)|,l < I < 2*m,0 < i < 2 2t+ V p )j 

< max {(2* - l)E,(2 t+1 - 1)X A } (35) 

where t = 0, 1, 2, .... Hence, if (EU holds, we have X A > X A . ■ 
Proposition 4.3 (A necessary condition for X A ): Let A Co (0) = A;i(0) = E, where (co,ci) is a companion pair. 
An achievable column correlation constraint Xf of M 2 ^ 2 , nW must satisfy 

X A > (2* - l)E. (36) 

Proof: (E3 implies that \A u ( t) (m)\ = (2* - 1)E for k = 2 2t+1 n^ - 1. Hence, Xf > max M {|A u w (Z)|, 1 < I < 
2*m, < i < 2 2 * +1 n(P)} > (2* - 1)E. ■ 
Corollary 4.2: Let ^4 Co (0) = -<4ci(0) = E and A^ = maxjA^jA^}, where (cq,ci) is a companion pair. For 



M 2'm,2*n( P >' * ^ l > when 



A o < (37) 



the minimum column correlation constraint 



Xi in (t,m) = (2 t -1)E (38) 

is achievable. 

Proof: If i31l holds, based on Proposition 4.2, M 2 '^ 2 * w satisfies the column correlation constraint A A = 
(2* - 1)£. On the other hand, Proposition 4.3 states that X A > (2* - l)E must hold. ■ 

Example 4.1: To construct the complex- valued complementary set matrices C( p ) and D^ p ^ with a column 
correlation constraint X A = 1, we can choose a companion pair Co = and ci = /*(co) = (j, — , j , +), 

where + denotes 1, — denotes —l,j = a/— 1 denotes the imaginary unit and j denotes — j, which satisfy A A < 1, 



i = 0, 1. Then, the companion matrix is 
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By employing length-extension JT31 . the complementary set matrix C^ 1 ) and its mate can be obtained as 
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Based on Corollary 4.1, and are complementary set matrices with a column correlation constraint X A = 1. 

Example 4.2: Let again Co = (+, j, and ci = /*(co) = (j, +), then, any sequence c in IZ^a satisfies 
\ A < A" 4 = 4. Starting with M| 4 = [C^ 1 ) D*^ 1 )] constructed in Example 4.1 and applying the size-extension d!7l i. 
we obtain M| g as shown in ( 1591 . Based on Theorem 4.1, M| g satisfies a column correlation constraint A" 4 = 4. 
On the other hand, Corollary 4.2 implies that M| 8 achieves its lower bound Xf = E = 4, since A^ = 1 < 

Example 4.3: Let us consider how to construct M| 8 satisfying a column correlation constraint S A = 12. Since 
t = 1 and m = 4, based on Proposition 4.1, a sufficient condition for = 12 is S A < 2, for E = 4. The 
companion pair (co,ci) in Examples 4.1-4.2 satisfies S A < Sq = 2, for i = 0,1. Thus, Mg 8 in must also 
satisfy a column correlation constraint S A = 12. Let {u^ 1 , < % < 31} denote column sequences of Mg 8 in d39l l. 
we can verify that max < i < 3l| = 12. 



+3 - 3-3-3 
3-3-3+3- 
-3 + 3+3+3 
3+3+3- 3 + 
-3+3+1 + 3 
1 + 3+ 3- 3 + 
+3 ~ 3 -1 ~ 3 
l-j-3+3- 



+3 - 3-3-3 
3-3-3+3- 
-3 +3+3+3 
3+3+3-3+ 
+1-1- 3-1 
3-1-1 + 1- 
-1 + 1 + 3+1 
3+1 + 1-1+ 



-3 + 3+ 3+3 
3+3+3 ~ 3 + 
+3 - 3- 3-3 
3-3-3+ 3- 
+1-1- 3-1 
3-1-1 + 1- 
-1 + 1 + 3+3 
3+1 + 1-1+ 



+3 - 3- 3- 3 
3-3 -3+3- 
-3 + 3+ 3+ 3 
3+3+3- 3 + 
+1~1 -3 -1 
3-1-1 + 1- 
-1 + 1 + 3+1 
3+3+1-1+ 



Remark: For the case t = 0, Corollary 4.1 implies that the correlation constraint for the companion pair is also 
the column correlation constraint of M 2 nn(p) = [C^ D( p )]. For t > 1, based on Proposition 4.1, small Sq may 
also help in reducing the column correlation constraint Sf. However, Corollary 4.2 implies that it is not necessary 
to search for the companion pair with smaller A^, once the lower bound A^ ifl (i, m) = (2* — 1)E has been achieved. 



B. Golay column sequences 

Theorem 4.2: Column sequences of complementary set matrices and D( p ) are Golay sequences, if and only 
if the companion sequences c and Ci are both Golay. 

Proof: Lemma 3.1 states that column sequences of and are either ±co or ±ci. Note that, a negation 
of a Golay sequence is also a Golay sequence. ■ 

We present a constructive method to obtain the companion pair from which an m by n complementary set matrix 
with Golay column sequences can be constructed, where m = 2 q+1 , n = 2 P+1 , p, q = 0, 1, 2.... 

Theorem 4.3: Let 



H 



(?) 

















2x2? 



h (5-l) h (9-l) 
n i,0 n i,l 



a i,l V n j,0 ) 



(42) 



where h-^ and are two row sequences of H^, i = 0, 1. The initial matrices are 



H 



(o) 



+ + 
+ - 



H 



(o) 



2x2 



+ - 
+ + 



2x2 



(43) 



Then, {h^ho^} and {h^,h^} are respectively Golay complementary pairs and, furthermore, /i(h^o) = 
and/ i (hj}) = -hg,forg = 0,l,2.... 

Proof: See Appendix C. ■ 

Example 4.4: Let q = 2, then Co = = (+ H 1 h) and Ci = h^ 2 Q = (H h H h ++). Based on 

Theorem4.3, {co,h ^} and {ci, ^} are, respectively, Golay complementary pairs, where hp { = (H 1 ) 

and = (+ + H 1 — ). Thus, the companion sequences cq and ci are Golay sequences. The length- 



extension ( fToT ) for p = 2 allows for constructing the following complementary set matrix 



c( 2 ) = 



+ 








+ 


— 




— 


— 


+ 


+ 


+ 


+ 










+ 


+ 


+ 




+ 


+ 


+ 




+ 


+ 


+ 


+ 









whose column sequences are Golay. Hence, the PAPR 



+ 


— 


+ 


+ 




+ 


— 


— 


+ 


— 


+ 


+ 


+ 




+ 


+ 




+ 






+ 




+ 


+ 


+ 




+ 


+ 


+ 




+ 


+ 



(44) 



J 8x8 



all column sequences of is at most two [17]. 



C. Number of zeros 



Proposition 4.4: Let the companion sequence Co be a length m sequence with z zeros, then, any column sequence 
of Mo! +1 „, ,„> contains 2 l z zeros. 

Proof: Based on Eqs. d2~Tl and d22t . Co and its companion ci have the same number of zeros. The number of 
zeros in each column sequence does not change after each length-extension operation. Each size-extension operation 
doubles the length of column sequences, as well as the number of zeros. ■ 

Example 4.5: In this example, we consider a ternary complementary set matrix and its mate with a column 

correlation constraint S A = 5. Let us set m = 8 and z = 1. We can find the companion pair Co = (H h++ +) 

and ci = /j(co) = ( h + H h 0) which satisfy S A < 5, i = 0, 1. The companion matrix is 



C <°) 



+ -- + + + 0+ 

--+++-+0 



(45) 



Using length-extension d!5l l. we extend C( ) as 



C (2) 



+ -- -- + -- 

-- + - + + + - 

- + + + + - + + 

+ + - + -- - + 

+ + - + -- - + 

+ -- -- + -- 

+ + 0- + 

+ 0- 0- 0- 



and D( 2 ) 



8x8 



- + + + - + -- 
+ + - + + + + - 
+ -- - + - + + 
-- + -- -- + 
-- + -- -- + 

- + + + - + -- 

0-0-0-0 + 

-0 + 0- 0- 



(46) 



8x8 



Based on Theorem 4.1, complementary set matrices and D^ 2 ) satisfy the column correlation constraint S A = 5. 
Furthermore, by setting t = in Proposition 4.4, we have that any column sequence of and D^ p ) has only 



one zero. 



V. Search for companion pairs 

A. Exhaustive search algorithm 

Let X(n,X) denote a subset of all sequences of X(n) which satisfy the correlation constraint A. For example, 
let B{2 t m) = {c | a € B, 1 < % < 2 t m}, then B(2 t m, S A ) = {c \ c £ B(2*m), S A < S A }, where B = {+1, -1}. 
Clearly, all column sequences of binary MO complementary set matrix 2 t n<p> w ^ tn a c °l umn correlation 

constraint S A can be found in B(2 t m, S A ). Hence, to construct a MO complementary set matrix M 2 „ t (p) 
whose column sequences are in X(n, A), we need to find a companion pair (cq,ci) such that T^c^ci ^ X(n,X) 
(see Lemma 3.2). 

Let us index all K = \X{m)\ sequences as Xj, for 1 < i < K. When a column correlation constraint A is given, 
desired companion pairs can be obtained by exhaustive computer search over X(m), as described in Table III. 

Note that the ACFs of sequences in TZ^} can be computed recursively using d3TT >. For binary sequences, we can 
simply check if x • y = to determine the companion pair. 

B. Minimum achievable column correlation constraint 

The exhaustive search algorithm in Table III can be easily modified to search for the companion pair with a 
minimum achievable column correlation constraint. However, the computing load is heavy, especially for large m 



TABLE III 

Exhaustive Search Algorithm 

j=0; 

for i = 1, 2, 3, K loop 
if iz£) C X(2 t m, A), 
3 = j + i; yj = x<; 
for I=i-l,j-2,...l. 

check if (yj,yO is a companion pair; 

end 

end 
end loops 



and t. Let t = in <l30b . the companion pair for the construction of M^ n(p) = [C( p ) Y)^] with a minimum 
achievable column correlation constraint is 

(Co, ci) = arg min {max {A x , A y } : (x, y) is a companion pair and x, y 6 X(m)} . (47) 

Based on Propositions 4.1-4.3, the above companion pair may also lead to the MO complementary set matrix 
2'nO) W ^ a re d uce d column correlation constraint for t > 1. Hence, in the following, we consider companion 
pairs leading to an achievable or a minimum achievable column correlation constraint A m j n (m) = A m j n (t = 0, m) 
for the case t = only. Table IV lists binary companion pairs with minimum achievable column correlation 
constraints A^ in (m) and S^ in (m). It can be observed that most of these companion pairs can achieve A^ in (m) 
and S^ in (m) simultaneously. 

When m is large, the exhaustive computer search is infeasible. Hence, the existence of a companion pair satisfying 
a correlation constraint is an important problem considered in the next section. 

VI. The existence of companion pairs 

The existence of a sequence of a desired correlation constraint has been studied in literature. For example, binary 
sequences with X A = 1 exist only for lengths 2, 3, 4, 5, 7, 11 and 13, and are called binary Barker sequences; binary 
m-sequences [23] with X p = 1 exist for length m = 2 l — 1, I = 2,3,4.... In this section, we exploit the correlation 
properties of companion pairs and analyze their existence for correlation constraints X A and \ p . 



TABLE IV 



Binary Companion Pairs 



m 






companion pair 


m 


\ A 




companion pair 


2 


l 


1 


- + 
++ 


12 


2 


8 


++++-+-++--+ 

-- + + + - + 


4 


l 


2 


+ 

- + ++ 


14 


2 


13 


+ - + + + -- + 

++--++-+-+ 


6 


2 


5 


- + + 

1 h+ 


16 


2 


/ 


- + + + + -- + -+ 

+--++-+-+ + 


8 


2 


6 


— 1 h+ 

+ + H 1 h 


16 


/ 


12 


+ H + + + + H 1 1 h 

- + + H h + + -H h-++ 


10 


2 


9 


+-+++++--+ 
+ + -- + -+ 


18 


2 


17 


+ + + + +-+--++--+ 

-+-++-+ + + --++ 



A. Correlation properties of the companion pair constructed using two arbitrary sequences 



In the following Cases 1 and 2 we illustrate how a companion pair of length m can be formed using two arbitrary 
sequences so and s\ of length m/2. We study the correlation properties of the companion pair (co, ci ) constructed 
from so and si. 

Case 1: Let Co = so <8> Si and ci = s* (8) (— Sq). Then ci is a companion of Co since ci = /*(co). 
In this case, ACFs of the companion pair can be expressed in terms of the ACFs of so and si and their 
crosscorrelation functions, 



A Co {l) = { 



Ac t (l) = { 



4s ,Si( 2 )+Aj ,Si( j 1 / 



I G odd 
Z G even 



-4 r*±i) 



A*, s ,( :: # i ), I e odd 



2 7 J1 s ,siV 2 

A So (^) + A S1 (|), I € even 



(48) 



(49) 



Pc (l) = { 



P SOiSl (^i) + P So , Sl (^^i), Zeodd 



P So (i) + P Si (i), 



I £ even 



P So (|)+P Sl (|), I € even 



(50) 



(51) 



where < I < m — 1. 



Lemma 6.1: Let c = s Si and = si <8> (— So), then 

A^<max{A^+A^,2A^ Si } 



and 



^ \r<max{\r+\r,2\^ si } 



Set - Ss + Sfh + ^,si 



where z = 0, 1. 

Proof: See Appendix C. 

Case 2: Let c = s S! and Ci = s*(— Sq). Then ci is a companion of c since ci = / c *(co). 
The aperiodic ACFs of the companion pair can be expressed as 

A So (l)+A Sl (l)+A So , Sl (l-%), 0<l<^ 

A So (l)+A Sl (l)-A So , Sl (^-l), 0<Z<§ 

— ^4s ,si (y — 

Lemma 6.2: Let c = Sq s i an d = Si(— s ), then 

Ac<A^+A^+A^ s , i = 0,l. 

Ci — »o 1 »i 1 So ,Si J * 

Proof: The proof is along the lines of the proof of Lemma 6.1. 



A ca {l)={ 



^<l<m 



(52) 



(53) 



(54) 



(55) 



(56) 



B. Existence 

Without loss of generality, we assume that are complex-valued sequences of length m/2 and ^4 S , (0) = P Si (0) = 
m/2, i = 0,1. 

Lemma 6.3 (Welch bound [25]): Let {sj,i = 0, 1, ...,K — 1}, denote a set of K complex-valued sequences of 
length N. If A Sz (0) = P Sz (0) = N for all i, then, 

/ K - 1 

(57) 
(58) 

V zi\ j\ — j\ — i 
where 



Pmax ^ -A^ 



if- 1 

NK-l 

K-l 
2NK -K-l 



max {Af , Af _ } 



An 



max {A^,A^ S .}. 



(59) 
(60) 



Proof: The proof can be found in [25]. 



1) Column correlation constraint X A : The following Theorems 6.1-6.2 restate the companion pair existence 
conditions from Theorems 4.1-4.2 in terms of the {s ,S!} pair existence conditions from Lemmas 6.1-6.2. 

Theorem 6.T. MO complementary set matrix M^ n(p) with a column correlation constraint X A exists if there 
exists a sequence pair {s ,S!} with A max = \X A . 

Proof: Theorem 4.1 states that MO complementary set matrix {p) satisfying the column correlation constraint 
X A exists, if and only if we can find a companion pair (cq,ci), such that, 

X A <X A , i = 0,l (61) 

Based on d52"l ). a sufficient condition for ( 16 It is, 

X A X A 
X A <—,i = 0,l and A s A o>Si < — (62) 



Based on d56b . 



< — , i = 0,l and A s A 0)Si < — (63) 

1 \A - 



Hence, by comparing d62l and d63l . we can set A max = ^A 

Proposition 6.1: Sequence pair {sq,si} of length y with A max = ^X A exists only if 



X A > 7 === (64) 

Proof: Let A max = \X A , K = 2, and N = ^ in (EU> of Lemma 6.3, then (f64j follows. ■ 
Corollary 6.1: Let Co = sq (g) si and ci = si ® (— Sq) be a binary companion pair of length m, and Uj denote 
column sequences of the constructed M 2 mn(p) , < i < 2n^ . Then, A^ < A^ < A^, where 



X A = max |a^«j , < i < 2n (p) } (65) 



A£ = max{A^ o + A^2A^ Si }, (66) 

= (67) 
V2m — 3 

Proof: Xyy is derived from Theorem 6.1 and d64b by noting that A^ is an integer for binary sequences. A^ 

follows from Lemma 6.1. ■ 
2) Column correlation constraint: X p 

Theorem 6.2: MO complementary set matrix (p) with a column correlation constraint A p exists, if there 
exists {so,si} of length m/2 with P ma x = ■ 

Proof: The proof follows along the lines of the proof of Theorem 6.1 and is omitted. ■ 



Proposition 6.2: A length m/2 sequence pair {sq^} with P max = \\ exists only if, 



A P > -j^= (68) 
\Jm — 1 

Proof: Setting P max = ±A P , K = 2 and N = § in <|57l> leads to ■ 

C. Achievable column correlation constraints 

Theorems 6.1-6.2 suggest searching for sequences so and si of length m/2 with good autocorrelation and 
crosscorrelation merits to form a companion pair with a small achievable column correlation constraint. Former is 
a long standing problem (e.g. see [18]— [21]). In [20], good binary sequence pairs with small \f o , and A^ s 
were found by using simulated annealing search algorithm, and were listed in Tables I and II. Based on Corollary 
6.1, we present A^y and A^ of their corresponding binary companion pairs in Table V, where the reference [20] 
indicates that data is obtained by using sequences from this reference. However, the cost function for the simulated 
annealing in [20] is not optimal in our case. We instead minimize the cost function 

/(s ,si) = max{A^ + A^,2A^ Si } (69) 

to obtain an improved Ag. 

In Table VI, sequence pairs {so,si} of length m/2 = 63,84 and 100 obtained using simulated annealing based 
on d69l are presented. The corresponding ACF merit A^ is calculated and compared to that of the sequence pairs 
from [20]. The proposed sequence pairs lead to companion pairs with an improved autocorrelation correlation merit. 



TABLE V 

Aw and Ag for Long Binary Seed Sequences 



m 


62 


74 


82 


106 


118 


122 


126 


134 


146 


158 


168 


182 


186 


200 


218 


240 




6 


7 


7 


8 


8 


8 


8 


9 


9 


9 


10 


10 


10 


11 


11 


11 


Ai [20] 


16 


18 


16 


18 


20 


22 


22 


22 


24 


24 


28 


24 


24 


28 


30 


28 


\ A 
A s 


13 


15 


15 


18 


18 


18 


19 


20 


22 


22 


24 


24 


24 


27 


28 


28 



VII. Conclusion 

We have considered a construction algorithm for MO complementary set matrices satisfying a column correlation 
constraint. The algorithm recursively constructs the MO complementary set matrix, starting from a companion pair. 



We relate correlation properties of column sequences to that of the companion pair and illustrate how to select an 
appropriate companion pair to satisfy a given column correlation constraint. We also reveal a method to construct the 
Golay companion pair which leads to the complementary set matrix with Golay column sequences. An exhaustive 
computer search algorithm is described which helps in searching for companion pairs with a minimum achievable 
column correlation constraint. Exhaustive search is infeasible for relatively long sequences. Hence, we instead 
suggest a strategy for finding companion pairs with a small, if not minimum, column correlation constraint. Based 
on properties of the companion pair, the strategy suggests a search for any two shorter sequences by minimizing a 
cost function in terms of their autocorrelation and crosscorrelation merits, from which the desired companion pair 
can be formed. An improved cost function is derived to further reduce the achievable column correlation constraint 
X A . By exploiting the well-known Welch bound, sufficient conditions for the existence of companion pairs are also 
derived for column correlation constraints X A and X p . 

We have left the general problem of finding MO complementary set matrices with a minimum column correlation 
constraint as an open question. An important step towards solving the general problem is to find new construction 
approaches for MO complementary set matrices. A design algorithm based on N-shift cross-orthogonal sequences 
can be found in [24]. However, their column correlation properties are intractable. 



TABLE VI 



Achievable X a for Long Binary Seed Sequences 



m 


merits 


s (or si ) 


si (or s ) 


126 


A s = 19 
Ail = 17 
As = 22 [20] 
Au = 17 [20] 


+-+-+--+-+-+++++++--+ 

--+-+++--+ + + + + + 

--++-++ + + - + + + + 


++-+++-++--++ +-++ 

-+++--++++--+++++-+-+ 
+ + + - + + + - + - + - + + 


168 


As = 24 
A„ = 20 
\i = 28 [20] 
Au = 21 [20] 


+--+--+++++ +-++-- 

+-++--++--++-++ ++ 

-+++-+-+ +++-++++-- 

+ - + + - + - + + + 


+ + + + + + + + - 

-- + + + - + + + + - + - + + 

-+-+--+++++ ++++++- 

++-++-+++-+-++-+--+-- 


200 


As = 27 
A;l = 23 
As = 28 [20] 
Ail = 25 [20] 


- + + + + + - + + + + + + -++ 

+ + + + - + - + - + + -- + -++ 

+ + + + + + + + - + - + -- + - + 

- + + + + + - + + + - + + --+- 

- + -- + + + + + - + - + - + ++ 


-- + + + - + + - + + + + --++ 

+ - + + -- + - + -- + -- + - + 

+ + -- + + + + -- + 

- + + + + + + + -- + - + + +- 

-++-+-++-++-+-+++--+ 



Appendix 

A. Proof of Lemma 2.1 

Let us first prove that {b*, —a*} is a mate of {a, b}. A proof for binary sequences can be found in Theorem 
11 of [3]. For complex-valued sequences, the complementarity of {b*, —a*} follows from 

Ap{l)+A_ p {l) = (A t (l))* + (A_^(l))* 

= A b (0+A_ a (0 
= A b (l) + A a (l) 
= 

for 1 < I < n — 1, where n denotes the sequence length. We further show that the pair {b*, —a*} is orthogonal to 
{a, b} in the complementary sense, as follows 



= (^.,t(0)*-\*(0 



= 



for every I. 

Refer to the proofs of Theorem 6 and Theorem 13 from [3]. If {ai,bi} is a complementary pair and {a2,b2} 
is one of its mates, then both {aia2, bib2} and {ai (g) a2, bi (g) b2} are complementary pairs. This completes the 
proof of Lemma 2.1. 

B. Proof of Theorem 4.3 



and 



H 



(o) 



H 



(o) 



" h (o) " 

n 0,0 




+ 


+ 


h (o) 

n 0,l 




+ 




" h (o) " 

n l,0 




+ 




h (o) 




+ 


+ 



2x2 



2x2 



It can be verified that {1iq q, hg ]} and {h^°Q, h^} are respectively Golay complementary pairs. Based on Lemma 
2.1, {h^QjhQ 9 ]} and {h^Q,h^}, q = 1,2,3..., constructed from d^2l are guaranteed to be Golay complementary 
pairs. 

We observe that /•(>$) = h$ and /i(hjg) = -hg. Let /^hg) = hg, /^hg) = -hg, then, 



/.(h^) = /,(h$ h$) = /tCbgg) (-/;(<{)) = h$ = . 
In a similar way, we have that /^(hQ 9 ^ 1 ^) = — h^i . This ends the proof. 

C. Proof of Lemma 6.1 

We give the proof for A^ < max { A^ + A£ , 2A^ Sj } and S£ < S^ Q + + S 1 ^ Si . Other proofs are similar. 
A£ = max{|P Co (/)|,l</<m-l} 

{772 777, ^ 

|p So (0 + p Sl (0l,i < i < - - 1; |P So , Sl (0 + Pso^H - 1)1.0 < i < ¥ - ij 

{777 777 ^ 

|P So (/)| + |P B1 (0I,1 < I < y " l;2|P So , Sl (0|,0 < Z < - - l) 

= max | Af 3 + A^ , 2A^ „ j 

I sq Si J So,Si J 



and 



m— 1 



5 c A o = (01 

z=i 

= ^ |^ SQjSl (o + 4, 0>8i h -1)1 + 2 1^(0 + ^(01 



1=0 1=1 

< Ei^wok E ^(oi + Ei^oWi + E 1^(01 

Z=0 Z=-2*+l 
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